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ABSTRACT 

One-dimensional models are formulated to describe pressure 
wave propagation through one- and two-component annular and mist 
flow patterns. The solutions include the effect of inertial interphase 
momentum in these configurations. In addition, the nature and impor- 
tance of interphase mass transfer is discussed for one-component 
systems . The models exhibit good agreement with the available data. 
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INTRODUCTION 

The phenomenon of pressure wave propagation in one- and two- 
component annular, annular-dispersed, and mist flow patterns has been 
the subject of several experimental investigations.^”^ The common ob- 
ject of these studies was to obtain pertinent information regarding the 
critical flow of such mixtures. The detailed results have shown conclu- 
sively that the propagation velocity and the critical flow rate cannot be 
directly related as they are for single phase systems. The data also 
show that homogeneous equilibrium models, such as that proposed by 
Steltz,^ are in error by as much as 200%. 

Fauske^ has shown that a realistic model for bubbly flows must in- 
clude the rate of interphase momentum transfer. Henryk presented a 
formulation which correlated the interphase momentum transfer in terms 
of the virtual mass associated with the discrete gaseous bubbles. This 
correlation exhibits good agreement with the available bubbly flow data 
in one- and two-component mixtures. In contrast, Henry et al.^) 
developed a propagation model for stratified and annular flows which as- 
sumes negligible rates for interphase heat, mass, and momentum trans- 
fer. The proposed solution is in excellent agreement with the two- 
component stratified data. The model was also compared to high quality, 
steam-water mist flows and agreed to within 10% of the available data, 
however, it did not predict the dependency upon the mixture quality which 
the experimental data exhibited. 

The one-dimensional analyses in this paper consider first the role 

of virtual mass in the interphase momentum transport. This concept is 
analyzed for both droplet and wavy- annular flows. Secondly, the inter- 
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phase mass transport rates for one-component mixtures are considered 
for rarefaction and compression waves. Based on geometrical and inter- 
face conditions, approximations are made to the thermodynamic equilib- 
rium relations to describe such processes. 

ANALYSIS 

The propagation of a pressure wave through a stationary gas- liquid 
mixture is a transient phenomenon when viewed in an Eulerian reference 
frame. As shown in Fig. 1, if the propagation is viewed from a Lagran- 
gian reference system moving with the wave, the wave front becomes 
stationary, and the previously stationary mixture will become a homoge- 
neous mixture (k = 1) which is approaching the wave front with the propa- 
gation velocity a Tp . As the amplitude of the wave becomes small, the 
gas and liquid velocities behind the wave become small such that 


This approximation enables one to express the one- 


dimensional continuity and momentum equations for the moving reference 


frame as: 


Hr -i du din 

a Tp — ap + (I - a)p ? + ap — £ + (1 - a)p ? = 0 (1) 

iF dZ L e L J g dZ 1 dZ 

and 

HP P H T i f du , du," 

- + a^ p — ap + (1- a)p 7 + 2a Tp ap — I + (1 - a)p 7 — - =0 (2) 

dZ iJ ^dZ L g £J Lr \_ g dZ L dZj 

It is assumed that a, p , and p, are of the form ffP(Z)] such that 

& 6 

_df_ _ _df dP 

dZ dP dZ 

then Eqs. (1) and (2) can be combined to give a general expression for the 
propagation velocity 
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The void fraction is given by 
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which can be differentiated to give 
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where the magnitude of the velocity ratio has been set equal to unity. 
Substituting Eq. (5) into Eq. (3) gives: 
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1 x(l - x) dP 

For the void fractions commensurate with annular, annular-dispersed, 
and mist flow patterns, the compressibility of the liquid phase can be 
neglected [(dp^/dP) = 0] and for the systems considered herein, the den- 
sity of the gaseous phase is negligible compared to that of the liquid, 

TWO-COMPONENT FLOWS 

In two-component systems such as air-water mixtures, the inter- 
phase mass transport is zero [(dx/dP) = 0]. This fact and the above ap- 
proximations simplify Eq. (6) to: 
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TP 


or + o/(l - a) 
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dp 


dP 


- + a( 1 - a)p 
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( 7 ) 
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In Ref. 10, Henry et al. presented a model for stratified and annular 
flows. As will be subsequently shown, this model, which assumed negli- 
gible rates of interphase heat, mass, and momentum transfer, is appli- 
cable to stratified and annular flows with smooth interfaces. This solu- 
tion will be outlined below since the derivation is pertinent to all the 
models presented herein. 

The formulation is one -dimensional, however, in a system such as 
an air-water mixture, the mismatch between the acoustical impedances 
of the phases is considerable and one must be aware of the two-dimensional 
aspects of the propagation. Davies states that, in annular flow, each 
phase will transmit a pressure disturbance at its own propagation velocity. 
Such behavior is shown in Fig. 2(a), and it is immediately apparent that 
this situation produces transverse pressure discontinuities which cannot 
exist. The relaxation of these discontinuities produces a transverse mo- 
tion similar to that experienced by pressure waves propagating through a 
liquid contained in a flexible pipe/ 12 ^ It is proposed that these transverse 
forces skew the initially planar pulse into a configuration like that shown 
in Fig. 2(b). When the depth of the liquid layer is much smaller than the 
length over which the pulse travels, the wave can be considered of con- 
stant shape and to travel with a single velocity. It is this velocity that the 
following derivations are concerned with, and all the approximations made 
must be considered in light of this two-dimensional flow pattern. 

Since the liquid is essentially incompressible, any temperature change 
due to the passage of a wave front is negligible. Therefore, any interphase 
heat transfer is dependent upon the gas behavior. Because of the large 
thermal capacity of the liquid phase (p ? c), the interphase heat transfer is 
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governed by transient conduction through the gas to the essentially con- 
stant temperature liquid. For the large gas volumes commensurate with 
typical stratified and annular flows, the thermal response of the gas can 
be neglected such that 


By definition 


^g = ^g 

dP yP 



u, 


( 8 ) 

(9) 


such that 


dP a Tp \dP dPy 


(10) 


For smooth stratified and annular systems, where each phase is continu- 
ous, the momentum equations for the separate phases can be expressed 
as 


dP 

dZ 
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and 
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— + P 7 a TP — = 0 
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Combining Eqs. (10), (11), and (12) gives 
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Equations (8) and (13) can be substituted into Eq. (7) to give: 
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The above ratio is essentially equal to unity for the mixtures of interest 
here (a > 0. 50). The prediction of Eq. (14) will hereafter be referred to 
as the "smooth annular” model. 

This separated flow model was used as an approximate solution for 
mist or droplet flows in Ref. 9. However, in such flows, the liquid phase 
is discrete and one must include the virtual mass of the discrete phase in 
the continuous medium. ^ If it is assumed that the droplet is spherical, 
the momentum equation for a given drop can be written as: 


dP 


dU , 


p u 

dZ L c dZ 


dZ 
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(15) 


If, in addition, it is assumed that the behavior of a given droplet is char- 
acteristic of all droplets in the flow field, Eq. (15) can be combined with 
the system momentum equation (Eq. (2)) and substituted into Eq. (10) to 
give: 


dk 

dP 



( 16 ) 


which for low pressures 


ip g «p l ) 


can be approximated by 
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(17) 


dk 

dP a TP p g (1 + 

The derivative dp /dP is dependent upon the rate of interphase heat 

o 

transfer. Here, as in the smooth annular model, it is assumed that, as 
a first approximation, the rate of heat conduction through the gas phase 
is negligible. 

The above information can be substituted into Eq. (7) to obtain an 
expression for pressure wave propagation in two-component droplet or 
mist flow configurations. 
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For a > 0.50, Eq. (18) is closely approximated by 


(19) 

a„ V 1 + a 

g v 

Most annular flows are characterized by a wavy interface. These 
surface waves also experience a virtual mass effect when there is rela- 
tive acceleration with respect to the gaseous core. The effect of surface 
waves will be investigated by analyzing a case where the wave amplitude 
is large compared to the minimum film thickness such as the configura- 
tion shown in Fig. 3. It is assumed that such a geometry can be approxi- 
mated by hemicylindrical filaments resting on the surface. For a cylin- 
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drical rod in cross flow, the coefficient of virtual mass is equal to unity. 
Therefore, the momentum equation for a given filament can be written as: 


dP 

dZ 


= Pl u l 


dZ 


g\ g 


dZ 



( 20 ) 


Following the procedure used for droplet flows, one can derive the follow- 
ing expression for the propagation velocity in the approximate flow con- 
figuration shown in Fig. 3. 


— = ( 21 ) 

a g 

This formulation will be referred to as the wuvy-annular model. 

COMPARISON WITH TWO-COMPONENT DATA 

Since each of the above models assumes that all the liquid is in the 
geometrical configuration considered in that model, the three solutions 
should bracket the results in the high void fraction region (a > 0. 50). 

Figure 4 compares the smooth annular model to the experimental 
data for quiescent stratified (or smooth interface) systems as reported 
in Ref. 10. (All two-component data presented herein has been cor- 
rected for 100% relative humidity. ) The agreement for these idealized 
conditions is excellent. 

Evans et al.^ and Hamilton^ reported experimental results for 
flows which are generally in the annular-dispersed or mist flow regimes. 
These data are shown in Fig. 5 and are in general agreement with the 
proposed models. Except at very low air velocities, films are relatively 
thin which minimizes the wave amplitude. Therefore, one would expect 
better agreement with the droplet model than with the wavy-annular ap- 
proach. The data and flow regime observations of Ref. I substantiate 
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such reasoning. The flow regime for the data point at a = 0. 57 should 
be approaching annular-slug transition and thus should contain large sur- 
face waves. Therefore, this particular point should be, and is, more 
closely approximated by the wavy annular model. 

The experimental data of Garrard^ is shown in Fig. 6(a). The void 
fractions reported in Ref. 3 are determined from film thickness measure- 
ments using wall mounted conductance probes. Such measurements are 

insensitive to liquid entrained in the gaseous core. On the basis of the 

('141 

measurements reported by Hewitt,' ' it seems likely that the vertically 
upward flows described in Ref. 3 would experience some entrainment'. 

To afford the reader a basis of comparison, the total liquid flow areas 
were estimated with the Lockhart-Martinelli^^ void fraction correla- 
tion. The data is replotted as a function of the calculated void fraction 
in Fig. 6(b). For the thin wall films reported by Garrard, one would ex- 
pect a relatively smooth annulus with entrainment in the core, thus, the 
data should be bracketed by the smooth annular and droplet models, which 
is indeed the case in Fig. 6(b). Since the annulus contributes no momen- 
tum transfer, the mist flow model and the data can be compared on the 
basis of the core void fraction. The void fraction in the dispersed core 
can be estimated by subtracting the liquid film area from the total liquid 

f 

area predicted by the Lockhart- Martinelli correlation. This comparison 
is given in Fig. 6(c) and it is seen that the experimental data are in good 
agreement with the mist flow prediction. 

In summary, propagation models have been presented for two- 
component smooth annular, wavy- annular, and mist flow patterns. The 
models emphasize the effect of the flow configuration on the rate of inter > 
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phase momentum transport, and all exhibit good agreement with experi- 
mental data from corresponding flow regimes . 

ONE-COMPONENT FLOWS 


The various momentum considerations discussed in the two- component 
section will be assumed valid for equivalent one-component flow patterns. 

A one-component mixture introduces an additional degree of freedom 
to the propagation phenomenon; namely, interphase mass transfer. As 
discussed in Ref. 9, the experimental data for bubbly flows indicates that 
the mass transfer rate in the wave front is negligible. However, this con- 
clusion may not be applicable to other flow regimes. 

Considerable effort has been directed toward analyzing pressure wave 
propagation in reacting gas mixtures/^’ ^ which is a problem closely 
allied with those considered herein. The major conclusion of these studies 
is that for any finite reaction rate, the wave front propagates in a frozen 
(no mass transfer, dx/dP = 0) manner. This particular conclusion is sub- 


stantiated by some very ingenious and detailed shock w 7 ave data taken in a 

( 18 ) 

supersonic nozzle. ' Another important conclusion of Refs. 16 and 17 
is that the amplitude of the wave front decays rapidly and the bulk of the 
wave travels at some velocity less than the frozen speed. It is shown in 


Ref. 17 that for times 


t » 


2X 



( 22 ) 


the bulk of the wave travels at a velocity close to that given by an equilib- 
rium mass transfer formulation. Therefore, an experimental investiga- 
tion which claims to have measured frontal velocities may actually be 



12 


measuring bulk wave velocities. 

When the frontal velocity is considered in light of the decaying ampli- 
tude and the two-dimensional aspects discussed above and illustrated in 
Fig. 2, it seems reasonable that the bulk velocity is the practical quantity 
to deal with when describing the overall one-dimensional behavior of a 
liquid-vapor system. Thus, the following solutions are intended to rep- 
resent the bulk wave behavior. 

As shown in Eq. (22), to determine the effects of mass transfer one 
must have a knowledge of the reaction (mass transfer) rate. To obtain 
such information for a two-phase system requires an understanding of 
the governing processes and the interfacial surface area. The latter 
requirement demands a detailed picture of the flow configuration which 
is the general unknown in two-phase flow. Therefore, one cannot present 
quantitative models, but by considering the nature of the processes and 
the general flow pattern, approximate solutions can be generated which 
give new insight into the overall one- dimensional behavior of one- 
component, two-phase systems. 

The nature of interphase mass transfer differs considerably between 
compression and rarefaction waves. A compression wave propagating 
through a mixture initially in equilibrium produces superheated vapor and 
subcooled liquid which are each individually stable states. To relax to 
equilibrium, the liquid must first be heated by the vapor which is a pro- 
cess essentially controlled by conduction in the vapor. Thus, the one- 
component compression wave case is very similar to two-component 
flows. In contrast, rarefaction waves generate superheated liquid and 
supersaturated vapor which are both metastable states and each may relax 
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independently of the other. Such a relaxation for the vapor requires con- 
densation which in turn requires a nucleation site. The liquid temperature 
is greater than that of the vapor, hence, the liquid phase cannot serve as 
the site. As illustrated by the condensation shocks witnessed in the flow 
of saturated vapor through converging-diverging nozzles, vapor can be 
considerably supersaturated before any spontaneous relaxation to equilib- 
rium is initiated. ^ Therefore, it will be assumed that in the wave front, 
no condensation exists and the vapor behaves isentropically as if it were 
superheated steam (y « 1.3). The superheated liquid vaporizes as it re- 
laxes to a stable equilibrium condition. Since there is a definite waiting 
time involved in the growth of vapor bubbles, it is assumed no bubbles 
are formed within the liquid phase. Thus, the vaporization occurs at the 
interface. Vaporization at the liquid interface yields a surface tempera- 
ture equal to the local saturation value which remains essentially constant 
regardless of the amount vaporized. The response of the liquid phase and 
the mass transfer is thus governed by the conduction response of the liquid 
volume with the imposed constant temperature boundary condition. 

As outlined above, the interphase mass transfer rate is governed by 
the transient conduction response of the vapor and the liquid for compres- 
sion and rarefaction waves, respectively. For identical geometries, the 
transient conduction response of water is more than two orders of magni- 
tude faster than that of steam. For typical annular and mist flows, 
the vapor volume fraction is one order of magnitude greater than that of 
the liquid. Hence, the transient response of the liquid is 1000 times faster 
than that of the vapor. Based on these crude estimates, it is assumed tha„t 
the transient conduction response of the vapor is negligible (dx/dP = 0). 
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Hence, the solutions for compression waves propagating through one- 
component annular, wavy-annular, and mist flows are identical to the two- 
component solutions presented above. (One possible exception will be 
discussed later. ) However, the mass transfer response in rarefaction 
waves requires further consideration. 

As was discussed above, the rate at which the mass transfer occurs 
is dependent on the nature of the process and the area available for 
transfer has been equated to the rate at which liquid conducts heat to the 
interface or 


q = -K Z A S 



( 23 ) 


where A is the interfacial surface area per pound mass of liquid and 
dT,/dy) g is the temperature gradient at the surface. For heat transfer 
purposes two geometries are considered herein - thick films (smooth 

and wave-annular flows) and small droplets (mist flows). The surface 

\ 

area for mist flows can be one to three orders of magnitude greater than 
that for the annular flows. Hence, if the interphase mass transfer is 
appreicable, the effect will be most noticeable in the mist flow regime. 

It is assumed, as a first approximation, that the mass transfer rate for 
films is negligible because of their comparatively small surface area and 
thick geometry, but that the rate of small droplets can be approximated 
by the equilibrium rate for vaporizing liquid. Therefore the rarefaction 
wave models for smooth annular and wavy-annular flows are identical to 
the two-component solutions. For mist flows, the total system and vapor 
entropies are assumed constant such that 


ds 


0 


d 


(1 - x)s f + xs 


0 


( 24 ) 
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In summary, considerations of the interphase mass transfer char- 
acteristics has lead to the following approximate solutions for the bulk 
wave velocity of pressure disturbances propagating through one- 
component mixtures. 

(1) The interphase mass transfer in compression waves is governed 
by the conduction response of the vapor w T hich is a comparatively slow 
process. Thus, it was assumed the mass transfer rate was negligible 
and the solutions for the three flow patterns considered are the same as 
the two-component models. 

(2) In rarefaction waves the liquid conduction response controls the 
mass transfer. This is a much faster process, however, it was assumed 
the small surface area and thick geometry of the annular flows severely 
restricted this response and the mass transfer rate was negligible. 

Hence, these two models are also the same as their two- component 
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COMPARISON WITH ONE-COMPONENT DATA 

The steam-water stratified data of Grolmes and Fauske^^ is com- 
pared to the smooth interface, no mass transfer model in Fig. 7. It is 
readily apparent that the compression waves exhibit excellent agreement 
with the analytical prediction. The rarefaction waves also are in good 
agreement with the model, however, there is a small but consistent dif- 
ference between the two. This small discrepancy is probably the result 
of some surface evaporation but since the conduction resistance of the 
liquid film is appreciable, the liquid mass transfer rate is considerably 
less than the equilibrium value and, for all practical purposes, negligible 
as was assumed. 

White and D’Arcy^^ measured the propagation velocities of rarefac- 
tion waves in steam- water annular flow systems. The mixer was designed 
to generate annular flows in a horizontal test section. No information is 
given about the shape of the interface in these flows, however, the smooth 
and wavy-annular models should bracket the data. As is illustrated in 
Fig. 8 this is generally the case. (The void fractions were calculated by 
the authors using a correlation given in Ref. 23. ) 

Figure 9 compares the experimental results for compression and 
rarefaction waves propagating through a single phase medium (air) and a 
50 percent quality steam-water mixture in mist flow. These data, taken 
by DeJong and Firey,^ clearly illustrate the difference in the propaga- 
tion phenomenon between the two media. The one-component mist flow 
models developed herein show excellent agreement with the measured 
velocities of small amplitude waves. 
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Figure 10 compares the models and the experimental data of England 
et al.^ for rarefaction waves propagating through steam-water mist 
flows. (Since the single phase propagation velocity is slightly dependent 
on the air content in the steam, which was not specified, a was taken 

o 

to be the average value measured in slightly superheated steam. ) Again 
there is generally good agreement between the model and the data. As 
stated in Ref. 5, at the lower qualities (x ^ 0.20), some liquid was form- 
ing on the wall, thus, decreasing the rate of interphase mass transfer. 
For such flow regimes, the data should be bracketed by the mist and an- 
nular predictions which is indeed the case. In fact, since the wall film 
contributes no mass transfer, one would expect that propagation veloci- 
ties for annular-dispersed flows could be calculated based on the mixture 
quality in the dispersed core. This is substantiated in Fig. 11 which ex- 
hibits the experimental results of Collingham and Firey^ for rarefaction 
waves propagating through annular-dispersed steam-water mixtures. 

The data are reported as functions of both gross and core qualities. As 
seen in Fig. 11(a), the data are bracketed by the mist and annular models 
as a function of gross quality, and, as shown in Fig. 11(b), the mist flow 
model affords a good prediction when based on the quality of the dispersed 
core. 

In Figs. 9(b), 10, and 11, the mist flow velocities for rarefaction 
waves are slightly less than the analytical predictions. As was discussed 
above, the bulk wave behavior is preceded by a low amplitude wave front 
which travels at the frozen velocity. The small interphase velocity dif- 
ferences set up behind this front could cause some viscous momentum 
transfer in the bulk behavior which would produce such lower velocities. 
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To summarize, considerations of the interphase mass transfer pro- 
cesses lead to the conclusion that such phenomenon are negligible in all 
cases except for rarefaction waves propagating through mist flows. A 
model which assumed negligible vapor condensation and equilibrium liquid 
vaporization was proposed for mist flows. 

One final point to be mentioned here is that the models and data dis- 
cussed above are characteristic of small pulses traveling through a mix- 
ture which is initially in equilibrium. Therefore, they are representative 
of neither sound waves nor critical flow and should not be construed as 
’’sonic" or "critical’’ velocities. As discussed in Ref. 9 for bubbly mix- 
tures and in Ref. 24 for mist flows, the velocity of sound in two-phase 
media is frequency dependent because of the relaxation times associated 
with the interphase transport of heat, mass, and momentum. It has also 
been shown in Refs. 25 and 26 that the velocities associated with one- 
component critical flows are characterized by considerable non equilibrium 
in the mass transfer process, and such flows cannot be treated by assum- 
ing negligible interphase mass transfer. 

COMMENTS ON THE DATA OF SEMENOV AND KOSTERIN 

Other than the annular flow results of White and D’Arcy, the only 

one-component data reported for the void fraction range 0. 50 to 0.90 are 

( 27 ) 

the steam-water compression wave results of Semenov and Kosterin. v ’ 

In light of the other experimental studies enumerated above, this investi- 
gation requires special attention. The experimental results, which were 
obtained in vertical-up flow, are shown in Fig. 12. For comparison, the 

air-water data given in Ref. 27 for a horizontal channel are also given. It 
is obvious that the air-water data does not agree with the results of Refs. 
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1-3, which raises questions regarding what flow configuration the data 
characterizes. No flow regime observations are given, but it is possible 
that some type of slug flow existed and the propagation was of the inertial 
type described in Refs. 10 and 28. For most of the void fraction range 
reported, this could also be true of one-component flows. However, for 
void fractions in the range from 0.90 to 1.0, a dispersed-annular con- 
figuration seems more likely than a slug pattern. Therefore, are the 
results characteristic of compression waves propagating through a one- 
component dispersed flow? 

As was discussed in the previous section, the vapor controlled mass 
transfer for compression waves is generally so small that it can be ne- 
glected. However, Eq. (22) shows the bulk wave behavior is also depen- 
dent upon the ratio of the frozen and equilibrium velocities. In the void 
fraction range investigated, this difference is quite large. Thus, perhaps 
the data of Ref. 27 satisfy Eq. (22) and, thus, are representative of the 
bulk of the wave traveling at the equilibrium velocity. The limiting mass 
transfer curves are shown in Fig. 12 (only virtual mass momentum trans- 
fer is considered) and the equilibrium mass transfer solutions describe 
the general character of the data. 

From the above discussion, it is apparent that a more definitive 
study of this regime is needed. Such a study should include not only 
rarefaction and compression wave data but also detailed flow regime 
observations. 

SUMMARY AND CONCLUSIONS 

One-dimensional models are developed for the propagation of pres- 
sure waves through one- and two-component mixtures in smooth annular, 
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wavy-annular, and mist flow configurations. The two-component models 
describe the nature of the interphase momentum transfer resulting from 
the virtual mass of the liquid configuration. In addition to momentum 
transfer, the one-component solutions incorporate assumptions based on 
the nature of the mass transfer processes in rarefaction and compression 
waves as well as the flow configuration. The various solutions presented 
exhibit generally good agreement with data characteristic of correspond- 
ing flow patterns. 

The following points can be concluded from this investigation, 

(1) If the phases are not configured in a smooth interface manner, 
the virtual mass of the liquid phase must be considered, 

(2) The vapor controlled mass transfer process for compression 
waves is considerably slower than the liquid controlled rarefaction pro- 
cess, and thus, the interphase mass transfer is more significant in 
rarefaction waves, 

(3) The interphase mass transfer is also strongly dependent upon 
the interfacial area which is determined by the flow pattern. Thus the 
interphase mass transfer is more significant in dispersed flows than in 
annular configurations, 

NOMENCLATURE 

A area 

a propagation velocity 
c specific heat 
f function 

K thermal conductivity 

k velocity ratio, u /u, 

& 6 
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P pressure 

q heat transferred per unit time per unit mass 

s entropy 

T temperature 

t time 

u velocity 

x quality 

y transverse length 

Z axial length 

a void fraction 

y isentropic exponent 

A. mass transfer relaxation time 

p density 

Subscripts: 

E equilibrium 

F frozen, dx/dP = 0 

f saturated liquid phase 

fg difference between liquid 

g gaseous phase 

l subcooled liquid phase 

o stagnation 

P constant pressure 

s surface area 

TP two-phase 
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Figure 1. - One-dimensional wave propagation velocity model. 
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Figure 3. - Approximate wave form employed to investigate 
momentum transfer in wavy-annular flow. 
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Figure 4. - Comparison of the smooth annular model with 
the stratified rarefaction wave data of reference 10. 
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Figure 5. - Comparison of the two-component models with 
experimental data of references 1 and 2. 
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Figure 6. - Comparison of models and data of reference 3. 
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Figure 7. - Comparison of the smooth interface model and 
the stratified data of reference 21. 
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Figure 8. - Comparison between the annular flow models 
' and the one-component rarefaction wave results of 
reference 22. 
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Figure 9. - Compression and rarefaction 
wave data . 
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Figure 10. - Comparison of the one-component models 
with the rarefaction wave results of reference 5. 
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Figure 12. - One- and two-component experimental results 
presented in reference 27. 



